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ABSTRACT 

We report the observation in the GeV band of the blazar 1ES 0229+200, which over recent years 
has become one of the primary sources used to put constraints on the Extragalactic Background Light 
(EBL) and Extragalactic Magnetic Field (EGMF). We derive constraints on both the EBL and EGMF 
from the combined Fcrmi-HESS data set taking into account the direct and cascade components of the 
source spectrum. We show that the limit on the EBL depends on the EGMF strength and vice versa. 
In particular, an EBL density twice as high as that derived by Franceschini et al. (2008) is allowed if 
the EGMF is strong enough. On the other hand, an EGMF strength as low as 6 x 10~ 18 G is allowed 
if the EBL density is at the level of the lower bound from the direct source counts. We present the 
combined EBL and EGMF limits as an exclusion plot in two dimensional parameter space: EGMF 
strength vs. EBL density. 

Subject headings: BL Lacertae objects: individual (1ES 0229+200) — cosmic background radiation 
— galaxies: magnetic fields 



1. INTRODUCTION 

Very-high-energy (VHE) 7-ray flux from distant 
blazars is absorbed on the way from the source to 
the Earth through its interaction wi th the Extragalac- 
tic Ba ckground Light (EBL) photons (|Gould fc Schrederl 
(|1967f )). The measurement of the induced distortions of 
the VHE 7-ray flux from distant hard-spectrum blazars 
by the effect of absorption on the EBL was used to derive 
const r aints on the EBL density (|Aharonian et all ([20061 
[2W^ : [QTret~aT1 (l20ll L 

The conventional derivation of the upper bound on 
the EBL from 7-ray observations adopts the assumption 
that the intrinsic powerlaw-type spectrum of the primary 
source (a distant blazar) is characterized by the slope 
dNj/dE ~ £+ r with T > 1.5. This assumption ap- 
pears reasonable in the framework of the most simple 
synchrotron-self-Compton (SSC) models for the broad 
band spectra of blazars. However, particular blazars 
considered for the derivations of EBL limits (the blazars 
with hardest intrinsic 7-ray spectra) may well not fit into 
this simplest SSC model framework, so that the assump- 
tion o f r > 1.5 might not be applicable fAharonian et al.l 
(120081) :lB5ttcher et all (120081) :lKatarzynskiet all (120061): 
iLefa et all (|2011al lbf): INeronov et al.l (|2011[ )). If the in- 
trinsic spectra of the blazars used for the derivation of 
the upper bound on the EBL density a re harder than 
r = 1.5, this upper limit is relaxed (Ma zin fe Raud 
(|2007l) ). 

Constraints on the intrinsic slope of the spectra 
of blazars can be obtained from the o bservations by 
the Fe rmi Large Area Telescope (LAT) (| At wood et al.l 
((2009T) ) in the energy band below ~ 100 GeV, where the 
effect of absorption on the EBL becomes negligible. How- 
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ever, the blazars used for the derivation of constraints on 
the EBL arc characterized by hard spectra, which makes 
it difficult to observe their flux below 100 GeV. In fact, 
the blazar 1ES 0229+2 00, which provides the tightest 
constraints on the EBL (jAharonian eTaTl (12001 ) is not 
listed in the catalo gue of sources det ected by LAT in 
two-year exposure (jAbdo et al.l (|2011D ). with only up- 
per li mits on the source flu x derived from the LAT 
data (INeronov fc Vovkl (|2010f ): iTavecchio et all ([2010D : 
Taylor et al. (120111 )) and a weak detection reported by 



Orr et all ([20"TT1) 



An additional difficulty for such constraints is that 
the spectrum of hard blazars might be composed of 
two contributions. Apart from the direct 7-ray emis- 
sion from the primary source, an additional contribu- 
tion is expected from the 7-ray cascade initiated in the 
inter galactic medium ( IGM) by the absorbed VHE 7 
rays (jAharonian et all (jl994f) ; IPlagal (|1995D ). The over- 
all flux and the spectral shape of the cascade contribu- 
tion are determined by t he strength of the extra g alactic 
magne t ic field (EGMF) (lAharonian et all (|2002D : IPlagal 
(fl995t ): [Neronov fc Semikozl (|2007l 120091 ). Uncertainty 
of the EGMF strength introduces an uncertainty of the 
importance of the cascade contribution and prevents the 
measurement of the slope of the intrinsic 7-ray spectrum 
of the source. In fact, the limits on the EBL derived up 
to now are based on an underlying assumption about the 
EGMF strength (the EGMF should be strong enough to 
suppress the cascade contribution up to the ^TeV energy 
band), which is not justified a-priori. 

If the assumption on the EGMF strength is relaxed, the 
7-ray data can be used to measure the EGMF strength. 
The electron-positron pairs, created as a result of the 
absorption of multi-TeV photons, up-scatter the Cosmic 
microwave background (CMB) as th ey cool, creating sec- 
ondar y emission in the GeV domain ([Neronov fe Semikozl 
(|2009D ). With the mean-free-path of TeV 7 rays be- 
ing D 1 ~ 100 TcV)- 1 Mpc, the production 
of these multi-TeV energy e + -e~ pairs occurs predom- 
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inantly in extragalactic space. Furthermore, if non- 
negligible magnetic fields arc present in the region these 
pairs are born into, they deflect significantly from their 
initial directions, resulting in the secondary GeV pho- 
tons having a reduced probability of reaching the ob- 
server. Thus, the apparent flux is suppressed at low en- 
ergies. Observational evidence for the presence of such 
a suppression places constraints on the E G MF strength 
to be > 10~ 17 G (iDermer et alj (|201~T1); iDolag et al.l 
cronov fc Vovkl (12010ft : iTavecchio et al.l (|2010t )7 
iTavlor et al.1 (|2011^ . 

This limit on the EGMF is derived assuming a cer- 
tai n level of EBL den s ity (a low EBL density derived 
by iFranceschini et ail (|2008fj is used in most of the 
publications on the EGMF in order to derive conser- 
vative bo unds). Howev e r, a v ariety of EBL model s 
do exist dGilmore et all (|2009f ): iKneiske et all (|2004l ): 
IStecker et al.l (|2006h ). which spans the full range of the 
present uncertainty in the EBL density. This uncertainty 
therefore introduces an uncertainty into the bounds on 
the EGMF derived from the 7-ray data. Higher (lower) 
EBL density leads to stronger(weaker) absorption of the 
primary 7-ray flux along with a stronger (weaker) cascade 
contribution to the observed flux in the GeV range. The 
suppression of such stronger(weaker) cascade contribu- 
tion would require a stronger (weaker) EGMF. 

Thus, the limits on the EBL derived from the 7-ray 
data depend on the assumptions made about the EGMF 
strength and vice versa. This implies that the correct 
procedure for the derivation of limits on the EBL from 
the 7-ray data should include marginalization over the 
possible EGMF values. Conversely, the correct proce- 
dure for the derivation of the limits on the EGMF should 
include marginalization over the possible EBL densities 
and spectra. 

Practical implementation of this correct procedure for 
the derivation of the EBL and EGMF bounds, however, 
present a challenging task since the quality of the 7-ray 
data is usually insufficient for the exploration of the en- 
tire EBL-EGMF parameter space. Exploration of this 
parameter space requires the measurement of the source 
spectra both in the TeV and GeV energy ranges. 

In this Letter we report the observations of the blazar 
1ES 0229+200 in the 1-300 GeV energy range using 
LAT's data with three years of exposure. The detec- 
tion of the source below 100 GeV provides the infor- 
mation necessary for the correct analysis of the EGMF- 
dependent upper bound on the EBL density and of the 
EBL-dcpcndent lower bound on the EGMF strength. We 
present such bounds in the form of a two-dimensional ex- 
clusion plot in the "EBL density" vs. "EGMF strength" 
parameter space. 

2. ANALYSIS 

In this work we use publicly available data of the LAT 
instrument collected over the period from August 2008 
till November 2011. We use the Pass 7 data and anal- 
yse them using the Fermi Science Tools v9r23pl software 
package, with the patches as of the November 1st 2011. 
We limit ourselves with the class 2 events, as is recom- 
mended by the Fermi/LAT tearr0. We select only pho- 

4 http: / /fermi. gsfc.nasa.gov/ssc/data/analysis/documcritation/ 
Cicerone/Cicerone_Data_Exploration/Data_prcparation.html 



tons with energies in the 1-300 GeV range for the analy- 
sis. During the spectral fitting, wc include in the analysis 
all the sourc e s liste d in the Fermi second year catalogue 
(jAbdo et al.l (]201 If )) within 5°-circle around the position 
of 1ES 0229+200. 

The fluxes in separate energy bins are computed with 
the spectral indices of all sources frozen at the values 
quoted in the Fermi/LAT second year catalogue; the 
spectral index of 1ES 0229+200 is fixed at the best- 
fit value in the 1-300 GeV band. For energies where 
1ES 0229+200 is not detected in the bin, we compute 
the 90%-confidence upper limit. 

To compute the "spectral butterfly" we scan the value 
of the likelihood in the "index-normalization" parame- 
ter space and select the corresponding 68% confidence 
region. 

We employ a Monte-Carlo description of the electro- 
magnetic cascades in order to determine the arriving 
spectra observe d for the cas e of d ifferent strength of 
the EGMF (see ITavlor et all (|201lD for more details). 
The shape of t he EBL spectrum is ass umed to follow 
that derived bv IFranceschini et al.1 (j2008l ). while the nor- 
malization of the EBL is left free. A lower bound on 
the normalization of the EBL from the dir e ct so urce 
counts (as summarized by iDominguez et all ([201 ID ) is 
~ 15% lower than the norm alization of the EBL model 
of IFranceschini et al.l (|2008f ). R ecent analysis of GeV to 
TeV spectra of several blazars (jOrr et al.l (|20 111 )), how- 
ever, suggested that the EBL level is somewhat higher, 
reachi ng ~ 60% of the EBL model of IFranceschini et al.1 
(|2008|) at the 2a confidence level. Both these bounds are 
taken into account in our analysis. 

In our calculations we consider the su ppression of t he 
cascade emission by time delay effects (jPlagal (|1995f) ). 
The TeV 7-ray emission from the source is observed to 
be stable on the time scale of > 3 yr, from the ini- 
tial H ESS observations of the source ()Aharonian et al.1 
(|2007f )) till the recent re-obs ervation by the Veritas tele- 
scope (|Dermer et al.l (|201lD ). We do not consider the 
alternative possibility for the suppression of the cascade 
flux due to the extended nature of the cascade emis- 
sion region. The generalization of our analysis to this 
case, however, is straightforward, the main difference be- 
ing that the suppression of the cascade emission due to 
the source extension would res ult in a somewhat h igher 
bound on the EGMF strength (|Tavlor et al.l (|2011l )). 

3. RESULTS 

1ES 0229+200 is detected by LAT in the 1-300 GeV 
energy band with significance ~ la. The test statistics 
(TS) value found in the likelihood analysis is TS = 45. 
Modelling the source spectrum as a powcrlaw we find the 
slope of the spectrum T = 1.36 ± 0.25 and normalization 
at 20 GeV (1.4 + 0.5) x 10" 15 (MeV-cm^sec)" 1 (at the 
68% confidence level). The spectrum of the source found 
from the LAT data is shown in Fig. Q] together with 
the HESS spectrum at higher energies. The source is 
not detected below ~ 3 GeV, only an upper limit on the 
source flux can be derived in this energy band. 

As is discussed above, the source spectrum in the 3- 
300 GeV energy band can have two contributions: the 
direct 7-ray signal from the primary source and emis- 
sion from the electromagnetic cascade developing in the 
IGM. It is not clear a-priori if the measured spectral 
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Fig. 1. — Top panel: the GeV-TeV SED under the assumption of 
a soft (r = 1.5) intrinsic source spectrum. Bottom panel: the GeV- 
TeV SED under the assumption of a hard (r = 1.2) intrinsic source 
spectrum. Lines for different values of EGMF strengths represent 
the sums of the intrinsic spectrum and the corresponding predicted 
cascade emission. In both cases the Fermi/LAT flux "butterfly" 
shown corresponds to the 68% confidence region. 

slope, consistent with T ~ 1.5, characterizes the intrin- 
sic source spectrum, the spectrum of the cascade com- 
ponent, or comprises a summed spectrum of the two 
(similar in strength) contributions. For instance, a spec- 
tral index hard e r then 1.5, as found in the analysis of 
iStecker fc SulTvl (|20"08| ) in the TeV band, would be in- 
dicative of a GeV spectrum which results from the sum 
of both the intrinsic spectrum and that of the cascade. 
Different possibilities for the dominance of one of the 
two components in the spectrum are illustrated in the 
two panels of Fig. [TJ In both models the normalization 
of the intrinsic spectrum is chosen to fit the HESS mea- 
surements in the TeV band. We also assume that the 
intrinsic source spectrum has a hi gh-energy cu t -off a t 
E cut = 5 TeV. As it was shown by iTavlor et al.l (|2011l ). 
this choice minimizes the strength of the cascade contri- 
bution in the Fermi/LAT energy band. 

In the upper panel the main contribution to the 3- 
300 GeV source flux is given by the direct flux of the 
primary source, shown by the thin solid line. This is pos- 




FlG. 2.— GeV-TeV SED under the assumption of a soft (r = 1.5) 
intrinsic spectrum for a different EBL scales. The EGMF strength 
of 10~ 16 G is adopted here, and the normalization of the spectra 
is chosen to fit the HESS data. The reference EBL scale is that of 
IFranceschini et al.l pOOST ). 



sible only if the cascade component is suppressed by the 
influence of a strong enough EGMF. If the EGMF is neg- 
ligible, the flux of the direct and cascade emission (thick 
solid line) will largely dominate over the direct emission. 
Strong EGMF (> 10~ 17 G) is needed to sufficiently sup- 
press the cascade emission down to the level of the error 
bars of the LAT measurements in the 3-300 GeV range. 

If the EGMF is weaker than - 3 x 10~ 17 G, the cas- 
cade emission provides the dominant contribution to the 
source spectrum, as is illustrated in the lower panel of 
Fig. [TJ The only possibility to make the LAT mea- 
surement consistent with observations is to assume that 
the intrinsic spectrum of the primary source has a slope 
harder than F = 1.5. The hardness of the intrinsic source 
spectrum depends on the EGMF strength. For the par- 
ticular example shown in the lower panel of Fig. [TJ the 
assumption that the EGMF strength B < 3 x 10~ 17 G 
imposes a constraint on the intrinsic source spectrum 
r < 1.2. In fact, if the intrinsic source spectrum is even 
harder, the intrinsic source flux contribution to the 3- 
300 GeV band flux becomes negligible and the flux is 
completely dominated by the cascade emission. 

The overall normalization of the cascade emission is 
determined by the density of the EBL. An increase of 
the EBL density leads to the stronger absorption of 
multi-TeV 7 rays and, consequently, to stronger cascade 
emission. To the contrary, reducing the EBL normal- 
ization down to the level of the lower bound from the 
direct source counts opens up the possibility of a weaker 
EGMF, down to ~ 6x 10~ 18 G. The effect of changing the 
EBL normalization is illustrated in Fig. [5J In this figure 
a spectral slope of T = 1.5 and EGMF of 10" 16 G have 
been adopted, and three different (66 %, 100% and 150% 
in term s of the EBL level reported bv IFranceschini et alJ 
(|2008t )) levels of EBL have been used. 

The maximal normalization of the EBL which can still 
be consistent with the data depends on the strength of 
the EGMF. Too strong an EBL can result in a large over- 
prediction of the strength of the cascade emission, even 
after taking into account the suppression of this emission 
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Fig. 3. — Allowed normalizations of the Franccschini et al. (2008) 
EBL spectrum, as a function of the EGMF strength. Dark grey 
shaded re g ions a re restrictions from [Dominguez et al. ( 20 ll|) an d 
Orr et al. (2011). The single stroke (cross) hatched region is ex- 
cluded by the analysis presented here at the 95% (99.7%) confi- 
dence level. 

by the EGMF effects. Thus, the upper bound on the EBL 
derivable from the 7-ray observations of 1ES 0229+200 
is EGMF-dependent. 

In order to find this bound, we compute the allowed 
ranges of the EBL normalization for a set of EGMF 
strengths and intrinsic spectrum spectral indices. The 
range of the EGMF strengths scanned over lay in the 
range 3 x 10~ 19 G to 10~ 14 G, while the spectral in- 
dices were varied within the range from 1.5 to 0. We 
then find the best-fit set of values in this "EGMF-EBL- 
r" parameter space, and chose an appropriate confidence 
region. The projection of this region onto the "EGMF- 
EBL" plane is shown in Fig. [3] and comprises the not 
hatched part of the plot. The hatched part thus rep- 
resents the EBL-dependent bound on the EGMF (or, 
equivalently, the EGMF-dependent bound on the EBL). 

One can see from this figure that for EGMF strengths 
B ~ 10~ 17 G the upper b ound sits at the level derived by 
iFranceschini et all (|2008f >. If the EGMF is at the level of 
~ 10 -15 G, the allowed EBL n ormalization is by a factor 
of 2 higher than that of th e IFranceschini et ail (|2008f ) 
and iDominguez et al.l (|2011| ) models. 

4. CONCLUSION 

Following our investigations into the newly detected 
blazar 1ES 0229+200 with Fermi/LAT we find that, if 



the IFranceschini et ail (|2008[) EBL level is adopted, an 
EGMF strength of at least 10 -17 G is required in order 
for the inevitable GeV cascade spectral component to 
be consistent with observations. However, we find that, 
more generally, the EGMF lower bound is dependent on 
the EBL level adopted, as demonstrated in Fig. [3] 

One should keep in mind that this bound is derived 
for large (As > 1 Mpc) correlation lengths of the 

— 1/2 

EGMF, and scales approximately as \ B ' for As < 
1 Mp c. Our result l ies in agreement wi t h pre vious find- 
ings (jTavlor et al.l pOill ); iHuan et al.l (|20TlT) L Under 
the assumption that the source is stable on a time-scale 
much long er than 3 years, a l arger b ound on the EGMF i s 
obtained (iNeronoy fe Vovkl (I20T0I); iTavlor et al.l ([20+1 : 
iDolag et all (|2010l ): iTavecchio et al.l (|2010l )L 



We also find that an intrinsic blazar spectral index of 
T = 1.5 is able to sit in agreement with observations of 
1ES 0229+200 by Fermi/LAT (r = 1.36 ± 0.25). Thus 
our results can find agreement with the underlying as- 
su mptions behind previou s EBL constraint calculations 
bv lAharonian et al.1 (|2006| ). However, the remaining un- 
certainty in the origin of the blazar's spectral shape in 
the GeV domain leaves open the possibility for a harder 
intrinsic spectrum (r < 1.5). 

We note, that in p r eviou s works of lAharonian et al.l 
(|2006h and lOrr et al.l (|2011| ) the authors made the im- 
plicit assumption for the absence of the cascade contri- 
bution in the GeV-TeV domain, which is equivalent to 
the presence of the strong EGMF with B > 10~ 15 G. 
As can be seen from Fig. [3l for such a strong field our 
result is also comp atible with the level of EBL found 
bv lOrr et all (j201lh . For weaker fields, ~ 10~ 17 G, we 
are als o compatible with the findings of lDominguez et al.1 
(|2011l ). However, for such weak fields a proper account 
should be given of the GeV-TeV cascade component con- 
tribution. 

Future observations with ground-base Cherenkov tele- 
scopes, such as MAGIC Stereo, HESS 2 and CTA, 
will dramatically improve the measurements of blazar 
spectra below ~ 100 GeV, allowing much better con- 
straints on the EGMF and EBL. In the particular case 
of 1ES 0229+200, if the intrinsic spectrum is hard, the 
presented EGMF lower bound may be transformed to a 
measurement, helping to better understand the nature of 
the extragalactic magnetic field. 

This work is supported in part by the SCOPES project 
No. IZ73Z0_128040 of the Swiss National Science Foun- 
dation. 
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